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ABSTRACT 



A model is constructed of a repairable item inventory system where 
on-hand stocks are resupplied both by repair and procurement. The basic 
properties of the model are analyzed. A Poisson failure process is 
assumed, and steady state probabilities are found for operation of the 
model under <R,Q> and<S-l,S> inventory policies. Further, a compound 
Poisson failure process under an <S-1 ,S> inventory policy is also 
analyzed. Measures of supply performance, such as expected number of 
backorders per unit time, are derived for the system under the various 

pol icies . 
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TABLE OF SYMBOLS AND ABBREVIATIONS 



C-| number of items in repair, X(t), at t=0 

C 2 number of items on order, Z(t), at t=0 

C 3 net inventory, Y(t), at t=0 

D delay betv/een item failure and establishment of a purchase 

order for a replacement; a random variable. 

D,(t) number of items repaired and returned to stock in the interval 
( 0 ,tj; a random variable. 

Dp(t) number of items ordered and received in the interval ( 0 , t] ; a 
random variable. 

f^ probability density function of procurement lead time. 

f^ probability density function of repair time. 

fy probability density function of batch size. 

F. u r c b 2 b i 1 i b y distribution 'Function ot 

Fp probability distribution function of repair time. 

H(t) the number of items on hand (in stock minus backorders) at any 
time t; a random variable. 

I time required to inspect an item; a random variable. 

ICP Inventory Control Point 

I(t) inventory position, i.e., on hand + in repair + on order, at 

any time t; a random variable. 

L mean procurement lead time. 

L time between placement of an order and receipt of material 

into stock; a random variable. 

N-j(t) number of failures entering repair during an interval of length t. 

N ? (t) number of failures which are lost to the system during an interval 

of length t. 

I 

N 2 (t) number of failures placed on order during an interval of length t. 

p probability an item which has failed is repairable 
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p-j probability a failed item is returned to the repair depot. 

P2 probability an item which arrives at the repair depot can 

actually be repaired. 

q probability an item which has failed is lost to the system. 

Q the amount of each order. 

r probability a batch of failures is repairable. 

R reorder level. 

R a mean repair time. 

Rq time required to repair an item; a random variable. 

<R,Q> type of inventory policy whereby an order is placed for Q items 

whenever inventory position first reaches x £ R. 

<s,S> type of inventory policy whereby an order is placed for S-x 
items whenever inventory position first reaches x <_ s. 

T-| in-transit delay time, customer to inspection; a random variable. 

T? in-transit delay time, repair depot to stock depot; a random 

yari a hi p 4 

U number of failures per batch; a random variable. 

1L 

U. number of failures in the 1 batch; a random variable. 

V ( t ) number of batches of failures in an interval (s,s+t]. 

V-|(t) number of batches of repairable failures in an interval (s , s+t] . 

Vp(t) number of batches which are lost to the system in an interval 

( s , s+tj . 

W(t) number of items which fail in an interval (s,s+t]. 

X ( t ) number of items in repair at any time t. 

Y(t) net inventory (on hand and in repair) at any time t. 

Z(t) number of items on order at any time t. 

A rate of a Poisson process. 
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I. INTRODUCTION 



Although inventory theory was one of the first areas to which 
operations research analysts turned their attention, most of the work 
has been concerned with the control of consumable items, i.e. items 
that are consumed in use. Until fairly recently the control of 
repairable items did not excite a great deal of effort. However, it 
has now been realized that, at least in military applications, repair- 
able or recoverable items are a major concern of inventory managers, 
both materially and financially. Repairables often account for a 
major portion of the total dollar investment even though the number of 
items may be small compared to the number of consumable items. Further, 
because repairable items are generally complex assemblies, they often 
require special storage, handling, packaging, and shipment methods. 
y One of the major di fficulties of applying inven to ry theory to 
consumable items is the extreme variability of item parameters, such as 
demand, procurement lead time, sources of supply, etc. Further, the 
very fact of very large numbers of items complicates an already diffi- 
cult task. Repairable items, on the other hand, have had a tendency 
to minimize some of these problems. The total number of items is much 
smaller and it seems somewhat easier to determine accurate parameters. 

For instance, a large number of repairable items may be classified 
electronic sub-components. Often, these can be further grouped and 
then treated as a group. The sources of supply and the sources of repair 
are generally completely known and limited in number. Repair lead times 
and procurement lead times may be often determined quite accurately with 
limited variances. Demand is generally more stable in the sense that if 
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a given repairable item is a sub-component of a larger unit, it is 
generally easy to find out how many larger units are in use, and this 
at least puts an upper bound on demand. 

On the other hand, the control of inventories of repairable items 
is complicated by the fact of two "sources of supply", i.e., repair or 
procurement. The purpose of this thesis is to analyze a model of a 
repairable item inventory system under some common inventory policies 
assuming Poisson and Compound Poisson demand processes. The question 
of optimizing the system under the various policies is not taken up, 
since inventory managers and other military decision-makers have not 
yet agreed as to the appropriate optimizing criteria. 
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II. A MODEL OF A REPAIRABLE ITEM INVENTORY SYSTEM 



The model used in this paper is a fairly close approximation of the 
actual repairable item inventory system in use in the Navy today. Almost 
every item stocked by Naval vessels and installations is assigned for 
technical and management control to some Inventory Control Point (ICP). 
Based on economic and technical considerations, the ICP will designate 
certain items as "repairable." These items are carried in stock, at 
various stock points or depots scattered throughout the Navy Supply 
System. Although for some items the ICP does not know how much of an 
item a certain stock depot has, repairable items are almost always 
"transaction reporting" items. This means that all issues, receipts, 
losses, etc. of that item at the Stock Depots are reported daily to the 
ICP, wmcn closely manages the item at a central level. 

Since items are often designated as repairable due to high procurement 
cost or critically low system stocks, the ICP is vitally interested in 
regaining all items which have failed. To facilitate this return of 
failed items, ships or stations requesting issue of repairable items are 
required to certify that a like number of failed items "has been or will 
be turned in." However, due to the ship's operating position, distance 
from stock points or repair depots, human procrastination, etc. there may 
be some delay in failed items arriving at the repair shop. Some items 

\ 

may be lost in transit or may never be turned in. Those that do arrive 
at the repair depot are inspected and some are found not-fit-for-repair 
and are discarded for salvage value or cannabalized for usable sub-parts. 
Those that can be repaired are eventually returned to the stock depots 
for reissue. Meanwhile the ICP, observing some items being lost to the 
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system, may decide to procure seme new items to bring the total system 
stock back to some desired level. The system just described is shown 
in Figure 1 . 




Figure 1. Actual Repairable Item Inventory System 



It should be noted that in actuality there usually will be several stock 
depots and probably at least two repair depots. 

However, this system can be considerably simplified with little or no 
loss of generality if it is considered from the point of view of the ICP. 
Typically, an ICP will know only v/hen items are issued, when items are 
received into stock from repair or procurement, and when failures are 
turned in to the repair depots. From this information, the total number 
in system stock or in repair at any time can be calculated easily. The 
ICP controls the procurement of new items and thus knows the number on 



order at any time. From these remarks, it can be seen that the various 
delays in the system, as seen from the ICP, are merely parts of either 
delay in repair or delay while on order. Therefore, it will be assumed 
that items which are repairable undergo a "repair time" which is the time 
from issue of the replacement item out of stock to receipt of the corre- 
sponding repaired item back into stock. Repair time in Figure 1 is 
therefore + I + R q + T^. Likewise, "procurement lead time" is the 
time between issue and return to stock of a replacement for items which 
are lost or not repairable and thus must be purchased. 

The basic quality of failed items is whether they are repairable or 
must be replaced. Assume failed items are turned into the repair depot 
with probability p-j and of those an item is repairable with probability 
P 2 , then failed items end up in repair with probability p^p^ and end up 

in nv'nmv'nniQni- wi-rh Ki 1 i i‘\/ f 1 — _ n ^ ~ ^ _ n Pi 

• ‘ * 1 ' — • • f- * — • .-s. i- • ~ * \ • • * i r I \ • i / -r illy 

• 1 1 ' 

Now, to simplify the model, assume that T.| , 1^, I, and D are equal to 
zero, or equivalently that the ICP has complete, instantaneous knowledge 
of the location and condition of all items, and can develop repair times 
and procurement lead times as given above. 

Let mean repair time be denoted by R^. Mean procurement lead time 
will be denoted by L. Assume that R^ has a probability density function 
f^ and probability distribution function F^. Likewise L has density 
f^ and distribution F L . Let p^p 2 = p and 1 - p^p 2 = q. Then the system 
described above can now bo modeled as shown in Figure 2. Assume that 
failure of any item is independent of the failure or non-failure of other 
items. Further assume that the number of customer requisitions in any 
period of time is a Poisson random variable with mean rate X. The number 
of items requested per requisition may be one each, in which case the 
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Delay = L 
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Delay = R ft 




Failures 
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Figure 2. Model of the Repairable Item Inventory System 



process is a simple Poisson process. If more than one item can be 
requester, rne process is a compound Poisson process. The number of 
items in repair at any time t is X(t). The number of items on order at 
any time t is Z(t). The number on hand, which is the number in stock 
minus the number of backorders, at any time t is denoted H(t). The 
net inventory of the system, defined as the number in stock plus in 
repair minus backorders at any time t is Y(t) = X(t) + H(t). 
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III. ANALYSIS OF THE MODEL UNDER POISSON DEMAND 



A. GENERAL ANALYSIS OF THE MODEL 

The number of failures in any interval (s,s+t] is assumed to be 
Poisson distributed with mean At. From the assumptions in Section II, 
it is seen that every failure results in either an additional item under- 
going repair or an additional item on order. Assuming that the 
classification of failures are mutually independent and independent of 
the process, it has been shown [Parzen 1962] that the number of items 
entering repair in the interval (s,s+t] is Poisson with mean rate Ap. 
Likewise, the number of failures which will be lost to the system in the 
interval (s,s+t] is Poisson with mean rate Aq. Further, the conditional 
distribution of the number of items entering repair in an interval 

given that n failures have occureo in that interval is binomial 
with parameters n and p. 

Now it can be shown that the number of failures entering repair in a 
given interval is independent of the number of failures which are lost 
to the system. Let N-j(t) be the number of failures entering repair 
during an interval of length t and N 2 (t) be the number of failures which 
are lost to the system during the same interval. Then the joint proba- 
bility mass function for both types of failures is: 



P(N-,(t) = x, N 2 (t) = n-x) 
= P (N, (t) = x, N(t) = n) 



= P (N, (t) = x | N ( t ) = n) P(N(t) = n) 
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x! (n-x) ! 




„ P x (n) x e-P u . (l-p) n - x (Xt) n - x e- (1 -P >xt 



x! 



(n-x).* 



, (pAt)VP Xt (gU) n - x e-q Xt 



x! (n-x)! 



= P(N 1 (t)=x) P(N 2 (t)=n-x) 



Thus N-j(t) and N 2 (t) are independent. 



At this point, before further analysis can be done, the processes 
of repair and procurement must be given some sort of formal definition. 
The number of items in the repair depot can be viewed as the number of 
busy servers in a queue, specifically an infinite server queue. For 
items being repaired; the idea of an infinite server cjue'ie Is not 
completely satisfactory. Certainly any repair depot set up to handle a 
continuing stream of failures will probably have the ability to work on 
more than one item at a time. However, depending on the Poisson failure 
process and the repair time distribution, the probability of large 
numbers of items entering or being in repair at any time may not be 
insignificant. Nevertheless, the assumption will be made that the 
repair process can be considered an infinite server queue, i.e. the 
repair shop immediately starts work on all items as soon as they are 
received. 

The other specification necessary for analysis is the form of the 
reorder policy. The inventory systems considered in this paper are all 
continuous review systems, that is, all transactions are recorded imme- 
diately, and thus, the inventory manager knows the exact value of stock 
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on hand, on order, or in repair at any time. A common continuous review 
policy is the <R,Q> policy, which specifies that a constant amount Q be 
ordered whenever the inventory position reaches R. The inventory posi- 
tion is defined as on hand plus on order plus in repair minus backorders. 
Under this policy the procurement process may be thought of as a queue 
where the number of busy servers at a time t represents the number of 
orders of size Q that are outstanding at that time. Another frequently 
used policy is the <s,S> policy, meaning order up to S whenever the 
inventory position falls below s. In the case where s = S-l , an order 
is placed everytime a failure cannot be repaired. Under the <s,S> policy, 
the procurement process is also comparable to an infinite server queue 
where the number of busy servers at a time t is the number of orders 
outstanding at that time. When items are demanded one at a time, the 
< R.Q > policy is the s?.rn° a ^ an < ^*S > nnlicv. 

Now let X(t) be the number of items undergoing repair at time t and 
Z(t) be the number of items on order at time t. Then 

X(t) = N-j (t) + C 1 - D 1 (t) 

Z(t) = N 2 '(t) + C 2 - D 2 (t) 



where C-j is X(0), C 2 is Z(0), D-j(t) is the number of items repaired and 
returned to stock in the interval (0,t], D 2 (t) is the number of items 

I 

ordered and received in the interval (0 , t ] , and N 2 (t) is the number of 
items lost to the system which have been placed on order by the time t. 

I 

N 2 (t) is seen to be related to N 2 (t) by: 



N 2 (t) = Q 



N 2 ( t ) 



Q JJ 



where jjn]i _ is the largest integer less than n. If items are reordered 

l 

every time one is lost to the system, i.e. Q=1 , then N 2 (t) = N 2 (t). 
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Otherwise ^ (t) <_ N 2 (t). Now both X(t) and Z(t) can be likened to the 
number of busy servers in a queue. D^(t) depends only on N, (t), the 

I 

repair queue discipline, and repair time distribution. N^ (t) and D ? (t) 
depend only on N 2 (t ) , the order queue discipline and procurement lead 
time distribution. Assuming independence between the two queue opera- 
tions, it is seen that X(t) and Z(t) are independent. 

Let H(t) be the number of items in stock ready for issue at time t, 
where H(t) is negative if backorders exist. The net inventory, Y(t) = 

X ( t ) + H(t), is the number of items in stock plus in repair minus back- 
orders. It is important for further analysis, that the independence of 
X(t) and Y(t) be shown, since H(t) will be obtained from X(t) and Y(t). 
This independence has been shown by Richards [Ref. 5] and it is instruc- 
tive to repeat the proof here: 



• . . exsii: I ! 10 thp irsnnncr in whinh transitions 0 m^dp in the two 
stochastic processes, (X(t), t _> 0) and (Y(t), t^O). 

Whenever a type-1 failure [into repair] occurs, H(t) decreases by 
one, but X(t) also increases by one for a net change of zero in Y(t). 

On the other hand, H(t) decreases by one while X(t) remains constant when 
a type-2 failure [procurement] occurs; thus a type-2 failure causes a net 
decrease of one in Y(t) but does not affect X(t). 

The arrival of a procurement increases Y(t) by the amount Q, while 
the repair level is unaffected. Conversely, X(t) decreases while H(t) 
increases by the like amount whenever an item leaves repair. Thus, 
repairs decrease X(t) but leave Y(t) unchanged." 

Then it can be seen that 



X(t) = C 1 + N-j (t) - D-j (t) 
Y(t) = C 3 - N^(t) + D 2 (t) 
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where = Y(0). As previously shown D^(t) depends only on N^(t) and 

I 

not on N^Ct). Conversely, N^ (t) and D^t) are functions of N^t) but 
not N^(t). Thus X ( t ) is a function only of N^(t) while Y(t) is a function 
only of N^lt). Since N^(t) and N^(t) have been shown to be independent, 
then X(t) and Y(t) are independent. 

B. AN <R,Q> POLICY WITH CONSTANT PROCUREMENT LEAD TIME 

Consider first the operation of the model under an <R,Q> policy, 
where the inventory position I(t) = X(t) + H(t) + Z(t) = Y(t) + Z(t). 

An order is placed for Q units whenever I(t) falls to R. Assume that 
all delay times other than repair time and procurement lead time are 
zero. Further assume that procurement lead time is a constant, L. Then 
the "net demand" placed on the system stock, Y(t), i.e. the number of 
items which must be replaced through procurement is a Poisson process of 
rate q. Thus the model from this point of view is merely an inventory 
system operating under an <R,Q> policy experiencing Poisson demands at 
rate q. This model has been solved by Hadley and Wnitin [Ref. 2]. 
Translating their notation into that used in this paper, the following 
are given: 

P[Y(t)=y] = 1 [F( R+l -y;XqL) - F(R+Q+l-y;AqL)] -*> <y 1 R+l 

= q- [1 - F(R+Q+l-y;AqL)] R+l < y < R+Q 

00 i -r 

where F(a;r) = 1 - £ — 

i=a 

i.e, the complementary cumulative of a Poisson distribution of rate r. 

Now 

Y(t) = X ( t ) + H(t) 
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thus 



H(t) = Y(t) - X(t) 

P[HCt)=h] = P[Y (t ) -X Ct ) =h] 

oo 

= z P[Y(t)=i+h,X(t)=i] 
i=0 

Since it has been shown that Y(t) and X(t) are independent 

oo 

P[H(t)=h] = z P£Y(t)=i+h] P[X(t)=i] 
i=0 

Note than an order is placed for Q units whenever the inventory position 
falls to R. Thus I(t) varies between R + 1 and R + Q. But 
R + 1 <_ I (t) < R + Q 

implies 

R+ 1 <_ Z ( t ) + Y ( t ) <_ R + Q 

Z(t), the number of items on order, is always non-negative, i.e. Z(t) 0 
Thus 

Z(t) + Y(t) < R + Q 

impl ies 

Y(t) < R + Q - Z(t) < R + Q 

or Max Y(t) = R + Q when Z(t) = 0 

Therefore P£Y(t)=y] = 0 for y > R + Q 

R+Q-h 

means P£H(t)=h] = 2 P£Y(t)=i+h] P£X(t)=i] 

i=0 

The probability that Y(t) = y has already been developed. Now P£X(t)=x] 
is merely an application of Palm's queuing theorem which states that for 
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an infinite server queue, with Poisson arrivals of rate A and service 
times distributed F(t) with mean M, the steady state probabilities of x 
busy servers are Poisson with rate am, i.e. in this case, 

tApR A )V X P R A 

P[X(t)=x] = * 

x! 

Since mean repair time is taken to be R^, and of course 



Then 



where 



E[X(t)] = ApR A 



R+Q-h 

P£H(t)=h] = E P[Y(t)=i+h] P[X(t)=i] 



i=0 



R+Q-h 

E 

i=0 

R+Q-h 

E 

1=0 



p(a;r) 



1 £F(R+1 -i -h ;AqL) -F(R+Q+1 -i -h ;AqL )] p( i ;XpR ft ) for i+h < R+l 

£1 - F(R+Q+l-i-h;AqL)] p ( i *,XpR^) for R+l £ i+h <_ R+Q 

= (r)V r 

a! 



Simp! ifying 



R+Q-h 



-1 II 

P£H(t)=h] = jr l £1 - F (R-h-i ;AqL) - 1 + F (R+Q-h-i ;AqL)] • p( i ;XpR A ) 
g i=0 

l R+ Q" h 

= ~r l [F(R+Q-h-i;AqL) p(i;XpR.) - F(R-h-i;AqL) p(i;XpR fl )] 

Q i=0 

, fR+Q-h R+h 

= E p(i;XpR A ) F ( R+Q-h-i ;AqL) - p(i ;XpR ft ) F (R-h-i ;XqL) 






i=0 



i-0 



Since H£t) = Y(t) - X(t) it is easy to find the expected value of stock 
on hand: 
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E[H(t)] = E[Y(t)] - E[X(t)] 



Richards showed that the expected value of net inventory for an <R,Q> 
policy is R + - m where m is the mean number of failures lost to 

the system in a procurement lead time. Thus 

E[H(t)] = R + E ~ ApR ft 

= R + ^r- " ^(pR A + PL) 

Now inventory position fluctuates between R+l and R+Q, i.e., 

R+l < Z(t) + Y (t) < R+Q 

Assume the system starts with Y(0) = R + Q and thus Z(0) = 0. Then when 
Y(t) drops to R, Z(t) becomes Q. If this order does not come in by the 
time Y(t) = R - Q, then Z(t) becomes 2Q and so forth. !he table below 



shows this relationship: 

Y(t) Z(t) 

R + Q 0 

r 1 , h 

R Q 

R -*Q + 1 Q 

R - Q 2Q 

« • 

• « 

R - 2Q + 1 2Q 

• • 

• • 

R - (n-l)Q nQ 

• • 

• • 

• • 

R-(nQ) + 1 nQ 
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Then Z(t) = nQ if and only if R-nQ+1 <_ Y(t) £ R-(n-l)Q n=l,2,... 

P[Z(t)=nQ] = P[R-nQ+l < Y(t) < R-nQ+Q] 

= F Y ( t )( R - n Q + Q) - F Y ) (R-nQ) 

x 

where F vm (x) = z P[Y(t)=y] 

UJ y=0 

Note that if Q=1 , i.e. reorder every time an item is lost to the system, 
this reduces to 

P[Z(t)=n] = F Y(t) (R-n+l) - F y(t) (R-n) 

= P£Y(t)=R-n+l] 

= ^ U(n;AqL) - F(n+l;AqL)] 

/■v „i \H^-AqL 

^ i t /\\j L / L' 

n! 

That is, when Q=1 , the number on order at any time t, like the number in 
repair, is Poisson distributed as expected by Palm's Theorem. 

Mow I(t) = Y (t) + Z(t) 

E[I(t)] = E[Y(t)] + E[Z(t)] 

r + r + - AqL + E[Z(t)] 

or E[Z(t)J= AqL 

Expressions have now been developed for X(t), Y(t), Z(t) and H(t) and 
the system is thus completely described. One measure of effectiveness 
that is almost always calculated for inventory models is the expected 
number of backorders at an arbitrary point in time. It is, of course 
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00 



E b p[H(t)=-b] 
b=l 

where b = the number of backorders. Under previous definition the 
quantity on hand is considered a negative value when backorders exist, 
thus H(t) = -b denotes b backorders on the records. 

Another well-known measure of effectiveness is the expected number 
of backorders per unit time which is also given by the expected number 
of failures minus the expected number of fills. The expected number 
of failures is just A, the rate of the Poisson failure process. Now 
a fill can occur only if H(t), stock on hand, is greater than zero. 

H(t) > 0 =>Y(t) - X(t) > 0 
=>Y(t) > X(t) 

Now Y (t ) is never greater than R + Q. If Y(t) > X (t ) , but X ( t ) > 0 
(always), then Y(t) > u. So 

p[fillj = p[Y(t) > X(t)] for 0 < Y(t) < R+Q 
Therefore : 

R+Q 

p[fill] = E p[Y (t) > X(t)|Y(t)=y] p[Y(t)=y] 

y=i 

p[Y(t) > X (t ) | Y(t)=y] = p[X(t) < y] 

y-1 (ApR A ) X e" XpR A 

= E 

x=0 x! 

= F(y-l;XpR A ) 

Then 

R+Q 

ptfillj = Z F(y-l;ApR A ) p[Y(t)=y] 

y=i 

R R+Q 

= E p[Y(t)=y] F(y-l;ApRj + E p[Y(t)=y] F(y-l;ApR A ) 
y=l y=R+l 
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R 1 

= 2 l [F(R+Q-y,XqL) - F(R-y,XqL)] F(y-l;XpR.) 

v=l q A 



R+Q 



+ Z tt F(R+Q-y;XqL) F(y-l;XpR fl ) 



y=R+l 



i Tr+q R 

= Q A ^ F(R+Q-y;XqL) F(y-l;XpR A ) - Z^F(R-y;XqL) F(y-l;XpR A ) 



Then the expected number of backorders is X - X • p[fi11] 




R+Q 

Z F ( R+Q-y ; XqL ) F(y-l;XpR.) 

y =l M 



R 






Z F(R-y;XqL) F(y-l;XpR A ) V . 



C. AN <S-1,S> POLICY WITH STOCHASTIC REPAIR AND PROCUREMENT LEAD TIMES 



The <S-1,S> policy has often been suggested as especially appropriate 
for repairable items since the object of the policy is to maintain the 



• ~ 4 - 1 ^ 

• IJ « * i >um r* i 



~ J i , 






i.c. when an item fails and cannot be 



repaired another is immediately placed on order. Assume that there is 
no lag between a non-repairable failure and the placing of an order. 
Then for this model 



X (t) + H(t) + Zf.t) = S for all t > 0 

The number in repair at any time t, X(t), and the number on order at any 
time t, Z(t) are modeled as infinite server queues, with arrival rates 
Xp and Xq respectively, where q = 1-p. In this model both repair time 
and procurement lead times are random variables with means R A and L 
respectively. Thus, Palm's queuing theorem is applicable to both X(t) 
and Z(t) giving 

(XpRjV XpR A 

P[X(t)=x] = - x=0 ,1 ,2 , . . . 

x! 
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p[z(t)=z] = ( x q L ) Ze ~ Xql 2=0, 1,2,... 

z: 

i.e. both X(t) and Z(t) are Poisson distributed. 
Now note that since X(t) + H(t) + Z(t) = S 



Y(t) + Z(t) = S 



then p[Y(t)=y] = p[Z(t)=S-y] y an integer, 1 y <. S 

_ (XqL) S - y e"^ C|L 
(S-y); 

Also p[H(t)=h] = p[Y(t)-X(t)=h] h an integer, £ y <_ S 



S 

= 2 p[Y(t)=y,X(t)=y-h] 

y=h 

C 

= ^ p[Y(t)=y] p[X(t)=y-h] 

y=h 

S (XqL) S -V XqL (A P Rj y - h e- A P R A 
p[n(t)=h] - ^ fs^TT (y-h): 



Feeney and Sherbrooke [Ref. 4], using a slightly different model, 
with an <S-1,S> policy, suggested the following three measures of supply 
performance (redefined for the present model): 

(1) R(S), called the ready rate, is the steady state probability 

that S or less units are in resupply. For the present model "in resupply" 
means in repair plus on order. 

„ R(S) = p[X(t)+Z(t) < S] 

S 

= £ p[X(t)+Z(t)=w] since both X(t) and Z(t) are non-negative. 

w=0 
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s 

= E 
w=0 



w 

E p(X(t)=x,Z(t)=w-x) 
x=0 



4 



S w 

= E E p[X(t)=x] p[Z(t)=w-x] 

w=0 x=0 



S w (X P R fl ) X e- X P R A (XqL) w - x e-^ L 

R(S) = E E - 

w=0 x=0 x! (w-x).' 



(2) F(S), called fills, is the expected number of demands per time 

period that can be filled immediately from stock on hand. As in the 
last section this is just the failure rate times the probability of a 
fill , i .e. X p[f i 1 1 ] 



p[f i 11] = p[H(t) > 0] 

= p[X(t) < Y(t)] 



= E p[X(t) < y|Y(t)=y] p[Y(t)=y] 

y=i 



S 

= ^ p[X(t) < y] p[Y(t)=y] 

y=l 



S-y e -XqL 



= E F(y-1 ;XpR A ) -^OL) 
y=l A (S-y)! 



S 

F (S) = XE F(y-l;XpR A ) 

y=i 



(XqL) 



S-y^-XqL 

(s-y): 



Now, it is well known that for a Poisson demand distribution the <R,Q> 
and <s,S> policies are the same, with R + Q = S and Q = S-s. In this 
case, with an <S-1 ,S> policy R + Q= S, R = S - 1, and Q = 1. In the 
previous section it was shown that for the <R,Q> model, 
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i \ R+ Q R 

PEfill] = n] ^ F(R+Q-y;AqL) F(y-1 ;ApR ) - Z F(R-y;AqL) F(y-l;ApR.) > 

Vy=l A y=! A 



Substituting the values for R+Q, R, and Q: 

S S-l 

Ptfill] = ^ F(S-y;AqL) F(y-l;ApR.) - Z F(S-y-l;AqL) F(y-l;ApR.) 
y=i A y=l A 

s-l 

= 2 [F(S-y) - F ( S-y— 1 ) ;AqL] F(y-1 ;ApR ) + F(0 ;AqL ) F(y-l;ApR.) 

y=l A A 



S-l 

= Z p(S-y;AqL) F(y-l;ApR.) + F ( 0 ; AqL ) F(y-l;ApR.) 
y= 1 M A 



Now F(0;AqL) = p(0;AqL) = p(S-S;AqL). 

Thus 

S 

p[f i 1 1 ] a Z p( S-y; AqL) F(y-l;ApR.) 

y=l A 

which agrees with the expression given above for the <S-1 ,S> policy. 

The expected number of backorders in a unit time is, again, 

A - A • p[fill] 

S 

= A - A Z F(y-l;ApR„) p(S-y;AqL) 

y=i A 

(3) A third measure of supply performance considered is S(S), the 
number of "units in service," defined as the expected number of units 
in routine resupply at a random point in time. Routine resupply is 
further defined as the number of units in resupply minus the number of 
units backordered. Again "in resupply" is taken to mean in repair plus 
on order. We note that if X ( t ) + Z(t) = r _< S, then H(t) 0, i.e., 
there are no backorders, and the number of units in routine resupply is r. 
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On the other hand if X(t) + Z(t) = r ^ S, H(t) < 0, i.e., there are 
backorders, and the number of units in routine resupply is S. 



Thus 

i 

s(s) = : 

r ; 



r 



S(S) = : 

r 



; r P[X(t)+Z(t)=r] + E S P[X(t)+Z(t)=r] 
=0 r=S 



: r^E P[X(t)=w] p[Z(t)=r-w] l + E s4 



=0 / w-0 



r=S 



E P[X(t)=w] p[Z(t)=r-w] 
w=0 



r (ApR.)V XpR A (Xql) r - W e- XqL 

: r-f I • — — — : 

: 0 ) w=0 w! (r-w)! 



00 1 r (ApR») e 

+ S E J E - 

r=S jw=0 wl 



w Q -ApRA ( X qL) r ' w e" AqL 



(r-w)! 
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IV. ANALYSIS OF THE MODEL UNDER COMPOUND POISSON DEMAND 



The model in the last section assumed that failures were a Poisson 
process of rate X, i.e. the number of failures in a time interval of 
length T has a Poisson probability mass function with parameter XT. 

This implied failures occurred one at a time, and thus each requisition 
for replacement was for exactly one item. 

In practice, the Poisson model may be quite reasonable for major 
repairable units, i.e. items of high cost and/or large, complex structure. 
However, there are smaller not so expensive items that are also economi- 
cally repairable and that may experience failures in batches, i.e. more 
than one at a time. One good example might be the printed circuit "cards" 
found in many modern computers, radars, and other electronic equipment. 
These cards are generally easily removed and installed, and under condi- 
tions of peak load, current surges, etc., several may fail at the same 
time. Generally, the only damage may be to the resistors, capacitors, 
transistors, etc. which are soldered onto the printed circuit, while the 
circuit card itself is undamaged and repair is feasible. 

These types of items lead to the consideration of failure processes 
in which the number of failures in a given time period is a compound 
Poisson process (W(t), t>0] . This process can also be represented by 

V (t) 

W(t) = £ U n 
n=l 



where { V ( t) , t>0} is a Poisson process which is the number of customer 
requisitions received and { , n=l,2,...} is a family of independent 
identically distributed random variables where U n is the number of items 
requested or the n th requisition. It is necessary that { V(t) , t>0} and 
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U n are assumed to be independent and for this model that means that 
the number of items per requisition is independent of which requisition 
it is, which is an entirely reasonable assumption. It is noted that the 
simple Poisson process is just a special case of a compound Poisson 
process, where the probability {U n =l,for any n} = 1 . The compound 
Poisson is an extremely rich family of distributions. Various seemingly 
quite different processes can be described merely by the proper choice 
of a compounding distribution. 

Assume that W(t) the number of items which fail in a unit time 
interval is a compound Poisson process with rate X and that U n is 
distributed as a random variable U, with probability mass function f 
and characteristic function <f>^. Then it has been shown [Parzen, p. 130] 
that the characteristic function of W(t) is 



■ • A ‘ wy ' x : 



> j 



Referring to the repairable item inventory model of this paper, 

assume that al 1 failure items in a particular batch are either repairable 

or all are non-repairable and result in a procurement. Further, let the 

probability that all items in a batch are repairable be r. Then the 

t h 

number of items to enter repair from the n u batch, R n> is 
U n with probability r 
0 with probability 1-r 

That is prob [R = u] = prob [U n = u] r u=l,2,3,... 
prob [R = 0] = 1-r 



The characteristic function of R n is 

00 

♦R- - EIe i0Rn ] = 1 e i0j p[R =j] 
0 j=0 " 



oo 



= p[R =0] + E i<t>j rn . n 
P n J j=1 e TJ p[R p =j] 

co 

= 1-r + r £ e 1 ^ p[U =j] 
j=l n 

= 1-r + r^ 
u n 



<f> R = 1 - r + 



Then ( t ) , the number of items entering repair in a time interval (0,t], 



V (t) 

N, (t) = Z R 
1 n=l n 



and from above 



_ e Xt(l-r+r4»u-l) 
= e Xtr (4>u “1 ) 



which is just the same compound Poisson process as failures with rate A 
replaced by rate Ar. Likewise, ^(t), the number of items which are 
placed on order in an interval (0,t], can be shown to be the same 
compound Poisson process with rate A(l-r). 

Now, let the number of batches entering repair be V^(t) and the 
number of batches which must be replaced by procurements be V^( t), so 
that V(t) = V-(t) + V 2 (t). The number of batches in the compound 
Poisson model corresponds to the number of items in the simple Poisson 
model. In Section II it was shown that the number of items entering 
repair, N^(t), and the number of items lost to the system, N 2 ( t ) . were 



independent random variables. In the compound Poisson process, uhis 
implies that V^t) and V 2 (t) are independent. Let N^t) be the number 
of failures in an interval (s, s + t) which can be repaired, and iet 
N 2 (t) be the number of failures in the same interval which are lost to 

the system. Then 

V(t) 

Ni(t)= E U nX] 

n=l 

V(t) 

N (t) = Z U n x 2 

n= 1 

fl if the n th batch goes to repair 
where X-) otherwise 

fo if the n th batch goes to repair 
*2 }1 otherwise 



distributed random variables, and a batch is either repairable or lost 
to the system. Thus N^t) and N 2 (t) are sums of independent identically 
distributed random variables, and, furthermore, no U n which is in the 
sum for N-| (t) is also in the sum for N 2 (t) and vice versa. Therefore, 

N ] (t) and N 2 (t) are independent. Further, X(t), the number of items in 
repair at any time t is a function only of ^(t) and the repair queue 
parameters. Likewise, Z(t), the number of items on order at any time t 
is a function only of N 2 (t), the procurement policy, and the procurement 
queue parameters. Thus X(t) and Z(t) are still independent for the 
compound Poisson case, assuming a batch goes completely to repair or is 
completely lost to the system. 

Again defining H(t) to be the number of items on hand at time t, and 
calling Y(t) = X(t) + H(t) the net inventory at time t, the same analysis 
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as in Section II shows that X(t) and Y(t) are independent. Using the 
notation of that section 

X(t) = C 1 + N-j (t) - D 1 Ct) 

Y(t) = C 3 - N 2 *Ct) + D 2 (t) 

I 

where C-j = X(0) , = Y(0), and N 2 (t) is the number of items which have 

been ordered in the interval (0,t], It has been shown that N^(t) and N 2 (t) 
are independent and it is assumed that D-j(t) is a function only of N^(t), 
i.e., that the number of items repaired in an interval of length t is a 
function only of the number of items which enter the repair queue and the 

I 

queue operation itself, and that N 2 (t) and ( t ) are functions of N 2 (t). 

Thus X(t) and Y(t) are functions of independent random variables and are 
themselves independent random variables. 

Assume vital the inventory system is governed by an tS 1 ,S> policy 
whereby every failure that cannot be repaired is immediately placed on 
order. This, of course, may result in an order for more than one item 
under the assumption that an entire batch is either repairable or not. 
Assume further that all items which enter repair or which are placed on 
order at the same time have the same repair time or procurement lead 
time respectively. The effect of these assumptions is that items which 
fail in the same batch are either repairable and are repaired as a batch , 
or that those which are not repairable must be ordered and are received 
into stock as a batch . While these assumptions may seem unduly restric- 
tive, in actual practice they are quite reasonable. Certainly items 
ordered at the same time are generally received at the same time. Further, 
repair often consists of a series of operations to be performed. It is 
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not unusual that the first operation is completed on all items in a 
batch before the second operation is started, etc., and upon completion 
of repair the items are returned to stock as a batch. 

Let be mean repair time and L be the mean procurement lead time. 
Let f u (x) be the "compounding" probability mass function, that is 
f u (x) = P{U = x} . 

Now Feeney and Sherbrooke have shown that under the above assumptions 
given a compound Poisson arrival density of rate A into a queue, the 
steady state probability of x units in the queue is the same compound 
Poisson density with rate AT, where T is the mean service time in the 
queue. Since the rate of items entering repair is Ar: 

Prob {x units in repair at any time) 

= Prob{X(t) = x) 

v fyv*D 

= i f 1 " (x) 

i=0 i! u 

1 * 

where f is the i-fold convolution of f u< 

Likewise, with the rate "into" procurement being A(l-r), 

Prob{z units on order at any time} 

= Prob(Z(t) = z) 

„ \ rA(l-r)LlV A(1 - r)L f i*( , 
i=0 i! U 

Then assuming the <S-1,S> policy, where an item is ordered every time it 
is lost to the system, the inventory position, which is equal to on hand 
+ in repair + on order - backorders, is always equal to S. In terms of 
the model 

H(t) + X (t ) + Z(t) = S 
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or Y(t) + Z(t) = S 



Then Prob (Y(t)=y) 

= Prob(S-Z(t)=y) 
= Prob(Z(t)=S-y} 



s-y 

E 

i=0 



[AQ-r)LlV A(1 ~ r)L 

i! 




(S-y) 



Furthermore, since Y(t) = H(t) + X(t) 



HCt) = Y (t) - X(t) < H(t) < S 



Then Prob (H(t)=h} 

= Prob(Y(t)-X(t)=h) 

S-h 

= E Prob(Y(t)=y+h,X(t)=y} 

w-A 
v Vy 

S-h 

= E Prob(Y(t)=y+h} Prob(X(t)=y} 

y=o 



and since Y(t) and X(t) are independent 

s-y-h [Ml-r)L]V X(1 - r)L y (XrR.)V XrR A jt 

Z f 1 , (S-y-h) Z - f ly) 

i=0 i! u i=0 u 

The system is now described for any time t 0. Again, the measures of 
supply performance suggested by Feeney and Sherbrooke are ready rate, 
fills and units in service. 

Ready rate, R(S) is the steady state probability that S or less units 
are in repair or on order, i.e. X ( t ) + Z(t) £ S. 

R(S) = Prob [X(t) + Z(t) <_ S] 

S 

= E Prob [X(t) + Z(t) = w] 

vr=0 



S-h 
= E 

y=o 



Prob [X(t)=x,Z(t)=w-x] 



S 

= E 

w =0 x =0 



w 

V 



S w 
E E 
w =0 x =0 



S w 

E E Prob [X (t ) =x] Prob [Z(t)=w-x] 

v /=0 x =0 



x (*rR A ) f e-* rR A w-x [X (1 -r )L ] 1 e ( 1 - r 5 L ., 

E f . 1 (x) E f 1 



i =0 



i =0 



(w-x) 



i : 



Fills, F(S), is the expected number of demands per time period for an 
item that can be filled immediately from stock on hand. Feeney and 
Sherbrooke reasoned that, "There will always be one fill if S-l or less 
units are in resupply [in repair plus on order] since each customer 
places an order for at least one unit. There will be a second fill if 
S-2 or less units are in resupply and the customer order is for two or 



- — ; .«*. • j -I ** 



Therefore, F(S) is the expected number of customers per time period, 
i.e. X, times the expected number of fills per customer. As shown above, 
this can be represented as 

p(S) = A ^ (probability a customer's demand is greater than 0) • 
(probability X(t) + Z(t) £ S-l) 

+ (probability a customer's demand is greater than 1 ) • 
(probability X(t) + Z(t) 5 ^ S-2) 



+ (probability a customer's damand is greater than S-l) 
(probability X(t) + Z(t) £ 0.)^ 

x(l • R(S-l) + [1 — f u (1 )]R(S-2) + [l-f u (l)- f u (2)]R(S-3) + 
... + H-f u O)-f u (2)-...-f u (S-l)] R(0) T 
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This is the same formula given by Feeney and Sherbrooke, except the ready 
rate has been adjusted for the double "sources" of supply, i.e. repair 
and procurement. Then the expected number of backorders in a unit time 
is, with probability one, the expected number of demands less the expected 
number of fills, i.e. 

E[number of backorders] = AE[U] - F(S) 

The other measure of supply performance is units in service, S(S), which 
is the expected number of units in routine resupply at a random point in 
time. There are r items in routine resupply if X(t) + Z(t) = r S, and 
there are S items in routine resupply if X(t) + Z(t) > S. Thus 
S oo 

S(S) = E r Prob[X(t)+Z(t)=r] + E S Prob[X(t)+Z(t)=r] 
r=l r=S+l 



S r 

E r E Prob[X(r)=w] Prob|_Z(l)^r-wj 
r=l w=0 



+ E S E Prob[X(t)=w] Prob[Z(t)=r-w] 
r=S+l w=0 



S r 
E r E 
r=l w=0 



XrR.)’e XrR A .* r-w [X 1-r )L]’e X(1 r)L .* 

A _f ( W ) E -fU r -v. 

i=0 i! u i=0 i! 



W 

z 



oo p 

+ S E E 

r=S+l w=0 



'w (XrR,) 1 e" xrR A 

z — f , M 

i=0 i : u 



r-w [X(l-r)L] i e" X(1 " r)L . 



E 

i=0 



i : 



f ' (r-w) 
u x 



36 



V. CONCLUSIONS 



The results developed in this paper are relatively straightforward 
extensions of models and results found by Hadley and Whitten and by 
Feeney and Sherbrooke. The model is a fairly good approximation of 
reality. The assumption of Poisson or compound Poisson demand processes 
has found much favor in inventory theory. The steady state probability 
of the number of items in repair, in stock, or on order, have been 
found, under various inventory policies. From these probabilities, 
formulas for several measures of supply performance have been derived. 
Other measures can be developed in a similar manner. Thus, given the 
steady state system probabilities and the derived measures of effective- 
ness, inventory managers can no w attempt to optimize the system under 
various criteria. 

Further analysis to improve the model may also be possible; for 
instance, the probability of a fill given a delay time T in which to 
make the fill would be a useful result. It may also be possible to 
relax the assumption that under compound Poisson demand, failure batches 
either go entirely to repair, or are lost to the system. It would seem 
more natural that each item of each batch may go into repair or may be 
lost to the system. 
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